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Abstract. We investigate the Kelvin-Helmholtz instabil- 
ity of stratified jets. The internal component (core) is 
made of a relativistic gas moving with a relativistic bulk 
speed. The second component (sheath or envelope) flows 
between the core and external gas with a nonrelativis- 
tic speed. Such a two-component jet describes a variety 
of possible astrophysical jet configurations like e.g. (1) a 
relativistic electron-positron beam penetrating a classical 
electron-proton disc wind or (2) a beam-cocoon structure. 
We perform a linear stability analysis of such a configu- 
ration in the hydrodynamic, plane-parallel, vortex-sheet 
approximation. The obtained solutions of the dispersion 
relation show very apparent differences with respect to 
the single-jet solutions. 

Due to the reflection of sound waves at the boundary 
between sheet and external gas, the growth rate as a func- 
tion of wavenumber presents a specific oscillation pattern. 
Overdense sheets can slow down the growth rate and con- 
tribute to stabilize the configuration. Moreover, we obtain 
the result that even for relatively small sheet widths the 
properties of sheet start to dominate the jet dynamics. 
Such effects could have important astrophysical implica- 
tions, for instance on the origin of the dichotomy between 
radio-loud and radio-quiet objects. 
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1. Introduction 

Kelvin-Helmholtz instability due to shearing induced by 
large velocity gradients is well known to develop in extra- 
galactic jets. This has been shown by several analytical 
and numerical works such as the ones reviewed for in- 
stance by Birkinshaw, 1991. These studies often consider 
the jet as a single fluid with one bulk velocity and a sin- 
gle interface made by a sheared layer with the external 
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medium. However there exist both observational and the- 
oretical pieces of evidence supporting the fact that com- 
ponents with different bulk velocities are present inside 
jets themselves. Among the observational data showing 
the presence of different velocities inside jets we can men- 
tion (i) the detection in some compact sources of VLBI 
components with different apparent speed, (ii) the super- 
luminal effect which proves that some jets are highly rel- 
ativistic at VLBI scale while properties of the extended 
component rather suggest classical or mildly relativistic 
large scale jets, (iii) the presence in some sources of en- 
trained gas which emits optical lines from which one can 
infer velocities likely smaller than the main bulk velocity 
of the jet. The peculiar morphology of the wide angle tail 
radio sources (WATs) appears also very difficult to explain 
(Eilek et al, 1984; O'Donoghue et al, 1990 and 1993) ex- 
cept if one assumes that matter with two different speeds 
is present in the jet, namely one relativistic component 
radiating before the inner hot spot and one component 
with slow speed radiating after the inner hot spot (Leahy, 
1984). Recent data on M87 provide the first direct mea- 
surement of apparent bulk velocity at the kiloparsec scale 
and show a complex velocity pattern with the presence of 
quite different bulk flow velocities along the jet (Biretta et 
al, 1995). In 3C273, the radio emitting jet appears longer 
and wider than the optical one. Bahcall et al (1995) sug- 
gest that there are actually two superposed components, 
a fast-moving inner jet surrounded by a slow-moving "co- 
coon" . Besides that, a Fanaroff- Riley type II radio source 
can also be considered as a stratified jet with two differ- 
ent layers corresponding to the inner real jet and to the 
backflow forming the surrounding cocoon. 

A quite new way of looking at the radio data appears 
as well in favour of the existence of several components 
in radio sources. Starting with the use of color-color dia- 
grams (Katz-Stone et al, 1993), Rudnick et al, 1994 and 
Katz-Stone and Rudnick (1994) show that it is reason- 
able to assume a distribution of radiating particles which 
is not a power-law throughout the sources and can then 
partially isolate the different parameters which contribute 
to the synchrotron brightness, namely the magnetic fleld 
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B, the number of radiating particles Nt and some fidu- 
cial Lorcntz factor 70 that characterizes their energy dis- 
tribution. Adding some knowledge concerning the shape 
of the radiating particle distribution, they can produce 
frequency-independent maps, proportional to NtB, JqB 
and Nt/jq. This powerful method provides a completely 
changed view of the sources in which new features are dis- 
covered. For instance, in the eastern lobe of Cygnus A, 
they detect an edge-brightened channel girdled by rings, 
which corresponds to a real enhancement of the radiat- 
ing particle density, likely located around the counter- 
jet (Katz-Stone, Rudnick, 1994). Recent developments of 
their technics lead them to suggest that jets in both FRI 
and FRII sources may have coaxial sheaths a few times 
wider than the jets themselves. From their "tomogra- 
phy" analysis of radio data which combine maps at differ- 
ent frequencies, they detect two-component structures in 
3C449 and 1231-1-674, with flat spectrum jets surrounded 
by sheaths of steep spectrum emission with different po- 
larization properties (Rudnick, 1995). It is not yet known 
whether such sheaths emanate from the inner jets or have 
been directly ejected from the central engines. However 
their presence deeply emphasizes the likelihood of multi- 
component jet models. 

Indeed several models of jet formation and propaga- 
tion reach the concept of stratified jets. First Chan and 
Henriksen, 1980, mentioned the "multilevel" structure as 
a basic modification to add to their self-similar jet model 
in order to explain why jets appear to have nozzles on 
very different scales. Smith and Raine, 1985, investigate a 
two-level configuration and explore the possibility to pro- 
duce coUimated outflows by the interaction of a nuclear 
wind from the very inner region of an active galactic nu- 
cleus with a Compton-heating-induced wind from an ac- 
cretion disc. Baker et al, 1988, deal with another kind of 
two-component model while studying the radiative prop- 
erties of a relativistic particle beam injected into an ex- 
tragalactic jet and the possibility of collective emission 
of radio waves. A completely different approach by Melia 
and Konigl, 1989, studies the Compton-drag deceleration 
of ultrarelativistic nuclear jets. It predicts the existence of 
a transverse gradient in the asymptotic bulk Lorentz fac- 
tor distribution of the particles in the radio jets. Konigl 
and Kartje (1994) also consider highly stratified winds 
from accretion discs, with inner ionized gas and outer 
dusty neutral outflow. Sol, Pelletier, Asseo, 1989, propose 
a two-component model taking into account a relativistic 
electron-positron beam extracted from the funnel of an 
accretion disc and streaming through a classical electron- 
proton coUimated wind coming out from all parts of the 
disc. This model provides a simple explanation to the ve- 
locity dilemma if VLBI jets and superluminal motion are 
related to the relativistic beam while kiloparsec scale radio 
features are associated to the slower coUimated wind. Sta- 
bility of such a configuration relatively to the excitation of 
Langmuir, Alfven and whistler waves has been found for 



electron-positron beams propagating along strong enough 

magnetic field with small enough Lorentz factor (Sol et al, 
1989; Achatz et al, 1990; Pelletier and Sol, 1992; Achatz 
and Schlickeiser, 1992). The condition for stability against 
excitation of Langmuir waves requires a strong longitudi- 
nal magnetic field B such that the electron gyrofrequency 
is larger than the plasma frequency in the wind, namely 



B > 3.2x10" 



in CGS units where n„ is the wind den- 



sity. The constraint on the Lorentz factor 7 comes from 
the necessity to avoid excitation of Alfven waves which 
imposes 7 < — 43. Therefore it appears possible to 

quench the beam-plasma instability and to ensure a two- 
component configuration stable from the point of view of 
microphysics. However, the question of large scale fluid 
instability of such two-component jet is still open. It is in- 
vestigated in the present work. As a first approach to the 
study of the Kelvin-Helmholtz stability of stratified jets, 
we consider here a core-sheet structure. It is the simplest 
"multilevel" configuration which appears in some jet mod- 
els (Smith and Raine, 1985; Baker et al, 1988; Sol et al, 
1989) and provides an approximation to the general case 
of jets surrounded by cocoons or sheaths. 
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Fig. 1. Schema of a stratified jet with core, sheet and external 
components in the unperturbed state. 



2. The dispersion relation 

We describe the jet with the 2-dimensional slab geometry, 
which provides an adequate description of the more real- 
istic cylindrical case with respect to the Kelvin-Helmholtz 
instability (Ferrari et al. 1982). For a large range of phys- 
ical parameters, the solutions for both the slab and cylin- 
drical cases have similar behaviour, apart from slight nu- 
merical differences. This similarity associates the symmet- 
rical and antisymmetrical perturbations of the slab jet 
with the pinching and helical perturbations of the cylindri- 
cal jet respectively. Only high order fluting modes do not 
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have counterparts in slab jets (Ferrari et al, 1982; Hardee, 

Norman, 1988). 

Both spatial and temporal stability analysis of Kelvin- 
Helmholtz instability have been carried out in the liter- 
ature. Spatial and convcctcd temporal growth rates have 
been shown to be not equal but still similar in the two 
approaches (Hardee, 1986; Norman, Hardee, 1988). This 
paper investigates the temporal behavior of the instabil- 
ity. We consider a core-sheet jet made of three layers as 
shown in Fig. 1 and describe the transition layers at all 
interfaces in the vortex-sheet approximation. 

The 2-dimensional core-sheet jet initially flows along 
the X-axis on the x—z plane. The unperturbed equilibrium 
state is defined as follows. Rc corresponds to the radius of 
the inner core, Rs to the radius of the surrounding sheet. 

{Vos ; ) is the bulk velocity in the core (respectively 
the sheet and the external medium) , Oc the sound velocity 
in the core, pc (ps, Pe) and Pc {ps, Pe) the proper density 
and pressure in the core (respectively the sheet and the 
external medium). We also introduce the density contrasts 

= Poo/ Pos and = Pojpo^, where index o stands 
for unperturbed zero-order values. The relativistic Lorentz 
factor of the core assumed to be a relativistic gas, the 
Mach number of the core and the Mach number of sheet 
are respectively 



7 

Me 

M, 



1/2 



-1/2 



as 



(1) 

(2) 
(3) 



We assume that the gas of the sheet moves with a non- 

rclativistic bulk velocity, Vo^ ^ c, and is described by 
a non relativistic equation of state with adiabatic index 
Tg = 5/3, however mildly relativistic sheet could also be 
taken into account within an analogous way. The external 
gas is also represented by a nonrelativistic gas with adi- 
abatic index Fg = 5/3, with 14 = in the unperturbed 
state. 

We shall use the relativistic equation of hydrodynam- 
ics in the form applied by Ferrari et al. (1978) for the 
relativistic core 



dt 



-(Ve-V)ye 



VPc + ^%=0(4) 



dt 



= (5) 



and for the non-relativistic sheet gas and external 
medium. 



^^''^ + (V,,eV)V,,e + ^ = 



dt 



Ps,e 



(6) 



dps,e 

dt 



V(p3,eV,,e) = . 



(7) 



The rest frame density of the relativistic core gas and 
the enthalpy are respectively 



mcUc + 



Pc 



lc2' 



Wc = Pc + 



(8) 
(9) 



where rric and ric are are the rest mass and the number 
density of gas. The sound speed of relativistic gas is (Taub 
1948, Landau & Lifshitz, 1959) 



(10) 



In the ultrarelativistic limit the adiabatic index Fc = 4/3, 
then 



{Vc-lf/^c- 



x/3- 



(11) 



The Lorentz factor of the core expressed by dimensionless 
quantities is 



7=(l-/3lM,2)-V2, 



(12) 



where (il^ = al/c^ = 1/3. 

For the nonrelativistic sheet and external gases we con- 
sider equations of state 



Ps,ePs^e"'' = Pos,e Pol','" = Constant 
with Ts,e = 5/3, hence 

Ps,e ^ 



(13) 



The sound speeds of the sheet and external gases are 



(14) 



^ s^ePs 



1/2 



(15) 



and the density contrasts can be expressed by sound speed 

ratios 



Po^ 

Po, 

Fea2 • 



(16) 
(17) 



From now on wc shall use dimensionless quantities 
with sizes scaled to Rc and time scaled to Rc/cloc and 
introduce the relative pressure perturbations 



9c,s,e — Pc,s,e/Po 

where Po = Poc = Pos = Poe 



(18) 
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The equations describing the relativistic core gas can while combining linearized p^ ) with (|2^) and for the 

sheet and with (|2^) and (pO|) for the external medium gives 



be written as 



r 2 



dt 



Vg, + PlV,^^Q (19) 
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^ + ^^^^^ 



= 



(20) 



and the equations of motion for the sheet and the external 
medium as 



dt 

dgP 
dt 

dVe 
dt 

dgP 
dt 



'Vg, = 



{VeV)Ve + VsVe 5e 9e = 



(21) 

(22) 
(23) 



The linearization of (|T^ and ( pl| ) leads to 



I dt 



^'^^ I dt 



dx , 



dx , 



dz 







(25) 



1 



1 



dt 

1 g^^. 
VsVe dt 



dx 



dx 



= 



(27) 
(28) 
(29) 
(30) 



in the rest frame of the external medium, where prime 
denotes first order perturbed quantities. Combining (p5|), 
(Eq) and linearized (EG) gives 



dz-^ 



d_ 

dt 



d_ 

' dx 



5c = 



(31) 



d 



2„' 



dz'^ 
and 

dz'^ 



d^ 

dx"^ 



dx'^ 



1 



92 



Fc dt^ 



5e = . 



5s = (32) 



(33) 



Let us assume in the core and the sheet perturbations 
of the form (for z > 0) 

6'c.s = [StsF^si^) + SZsF-si^)^ expz(/c,a; - cut) + c.c. (34) 

where F^^{z) — exp(±i fcc^,s^z) allow to describe waves 
propagating in opposite z-directions. Perturbations in the 
external medium are of the form 



(5g = d'^F^{z) exp i{kxX — Lot) + c.c. 



(35) 



where F^{z) — exp i k^^z represents only outgoing waves 
(24) for positive z. Here is the parallel wavenumber longi- 
tudinal to the jet and k^^ (respectively k^^ and k,,^) is 
the transverse wavenumber perpendicular to the jet axis 
in the core region (respectively in the sheet and external 
regions). The transverse wavenumbers fcc^, /cs, and /cg, 
are complex since the perturbations we consider represent 
surface waves. 

For a temporal stability analysis, k^ is real and lo com- 
(26) plex. Equations (31), (El) and (M) lead to 



ke.. — 



Xo 



1/2 



1/2 



(36) 
(37) 

(38) 



where uJc^ — j{uj — Mck^) and k^^ = "/{kx — l^a^McUj) are 
the frequency and wavenumber in the rest frame of the 
core fluid and — uj — ^/vsT sMskx is the frequency in 
the rest frame of the sheet fluid. From equations (10), (27) 
and (p9) we deduce the relations 



(39) 
(40) 
(41) 













±iysks. 


G'± 


V'+ = 


^s^eke^ 
LU 





Symmetry conditions of the problem allow to consider 
only half-space with z > 0. The components (-I-) and 
(— ) propagate respectively in the positive and negative 
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z-direction. Our perturbation takes into account only out- 
going solution in the external medium. 

The dispersion relation of the system is obtained from 
the boundary conditions at the internal (core/sheet) and 
external (sheet/external gas) interfaces, namely pressure 
equilibrium and equality between the ^-component of the 
gas displacement and the transversal displacement of con- 
tact surface. In the reference frame of the core gas, this 
condition writes 



d 



h'-fy = 



(42) 



We call h'^ and h'^ the transversal displacements of 
core/sheet and sheet /external surfaces. The exponent (c) 
stands for values determined in the core gas rest frame. 
From Lorentz transformation 



1 - I31{M, + v')M, 



(43) 



which reduces to v, 
Moreover 



ic)' 



d 



d 



ox 



jv'^ in the linear approximation. 



(44) 



and h'^, the displacement of surface, is Lorentz invariant 
as well as the pressure g'^. Boundary conditions at the 
internal interface located at z = 1 become 



(45) 



d 



d 







(46) 

9c = g's (47) 
and at the external interface located at z = i? (= Rs/Rc), 



where we have used the fact that by definition 
F+{R)F-{1)/F-{R) = F+{2R- 1) and have imposed 
the conditions of symmetry or antisymmetry of the solu- 
tions, — eG+ with e ±1. We introduce the complex 
normal acoustic impedances (Payne & Cohn 1985) for the 
core, sheet and external medium, which after elimination 
of factors constant across all the components can be ex- 
pressed as 



Zs 
Ze 



(54) 
(55) 
(56) 



The reflection and transmission coefficients at core-sheet 
and sheet-external gas interfaces are respectively defined 
as 



7^c 

Us 



Zs ^ Zc 



Zs 

Zc 



Zc 
Zs 



Zc^Zs 



Tc, 
Ts,- 



2Zs 



Zs + Zc 

IZc 
Zc + Zs 



(57) 
(58) 



When there is no reflection at this external boundary 
(interface sheet-external), Tl-sc = and the effect of the 
external gas disappears. Equations (^ reduce to the case 
of a single jet with internal medium corresponding to the 
core component and external medium to the sheet com- 
ponent. Another way to cancel the effect of external gas 
is to assume a very large sheet since the factor of TZse in 
( ^ vanishes as exp[— 2i?Im(ks^)] for R >> 1. A strati- 
fied jet with infinite sheet therefore corresponds to a single 
jet with only core and sheet components. Conversely, for 
R = 1, internal and external interfaces coincide and equa- 
tions (HI) reduce again to the case of a single jet, but then 
with only core component and external gas. 



/ 

9e ■ 



(48) 

(49) 
(50) 



We can take the forms of the transversal displacements 
as follows 



h'c s — ^c,s exp i {kxX — ujt) + c.c. 

After substitutions and elimination of ^c,s, 
general dispersion relation for core-sheet jets 

Zs\ F+{1) + TZscF+{2R - 1) ^ Fc+{1) + ef '(1) 
Zc) Fs+{l)-TZscFs+{2R-l) F+{l)-eFc{l) 

which leads to 



(51) 

we reach the 



(52) 



Zc 



i + nsef^^2}!lll 



coth ikc^ for e = 1 

''or (53) 
th ikc^ for £ = — 1 



3. Numerical results and their physical interpreta- 
tion 



The method used to solve the dispersion relation ( |53D is 
based on the Newton- Raphson method (Press et al, 1992). 
Different sets of parameters have been investigated. They 
all lead to the same qualitative results. 

1. The presence of the sheet induces an oscillating pattern 
on the growth rate diagram 

2. The stability properties are dominated by the sheet 
even if its thickness is comparable with the core radius. 

3. The amplitude of oscillations of growth rate diminish 
with increasing sheet thickness. 

This is well illustrated in Figs. 2 to 4 which show the tem- 
poral growth rate of the antisymmetric fundamental, first 
and second reflection modes of core and sheet jets. 

The apparent oscillating pattern is a new feature 
caused by the reflections of acoustic waves at the inter- 
face between sheet and external gas. These reflections are 
represented by the reflection coefficient TZse in dispersion 
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Fig. 2. Temporal growtli rate of the antisymmetric solution for an underdense core and sheet jet, for fundamental mode (top), 
first reflection mode (middle) and second reflection mode (bottom). R — Rs/Rc ~ "2, 5 and 10 for first, second and third columns. 
Dotted lines correspond to a single jet case with no sheet and R — 1, while dashed lines correspond to a single jet case with no 
external component and R — oo (infinite sheet). Parameters are 7 — 10, A'ls = 0, Vc = pel Pa ~ 0.01 and = Ps/pe ~ 0.1 



relation (|53|). The wave radiated at the core-sheet (cs) in- 
terface toward the sheet layer is reflected from the sheet- 
external (se) interface and next returns to the (cs) inter- 
face as shown in Fig. 5. Thus, one can expect resonances 
due to the presence of the sheet component if 



A(j) = 2mT, 



(59) 



where is the change of phase of sound wave in sheet 
on the path between cs, se and cs interfaces and n is an 
integer number. Introducing the length of path of the wave 
between cs and se interfaces (see Fig. 5) 



L = 



R-1 



cos a 



(60) 



where a is the wave propagation angle defined by tga 
kx/ks^, and the wavelength 



A., 



kl + fc2^ 



(61) 



we can write the resonant condition (159) in explicit form 



2L — riA, 



(62) 



The above interpretation of the oscillating pattern is 
demonstrated in Fig. 6., where we presented again the so- 
lution of Fig. 4 in the linear scale. The superimposed high 
bars represent these values of k^ (one should remember 
that ks^ is a function of kx) for which the condition (62) 
is fulfilled. For comparison, the small bars represent these 
kr for which we have an antiresonance 



2L 



(63) 



We can easily notice a coincidence of the maxima of 
growth rate with high bars and minima with small bars, 
what proves the validity of our physical interpretation. In 
a case of external gas which is significantly lighter then 
the sheet gas we notice TZse — —1 (the change of phase by 
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Fig. 3. Same as Fig. 2 but here Vc = 0.01 and Vs 
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7r) resulting in a reversed relation between bars and the 
maxima of growth rate. This effect is remarkable compar- 
ing Figs. 2 and 3, where maxima and minima in both the 
figures are inverted. Figure 5 shows a clear coincidence of 
maxima, minima and bars described above. Varying the 
parameters of the core sheet configuration we observe how- 
ever that the bars can be shifted uniformly with respect to 
maxima and minima toward smaller or larger values. 
This is not a surprise since all the reflection and trans- 
mission coefHcients are complex quantities in general and 
the additional phase shifts are natural for each reflection 
at the interfaces. 

The entire interaction of sound waves at the core-sheet 
interface appears to be very complicated. The normal 
acoustic impedances Zc and Zg depend on both real and 
oscillating imaginary parts of (jJ and fcc,s^. As Payne & 
Cohn (1985) shown the instability coincides with vanish- 
ing (in practice small) values of denominator in expres- 
sion for the reflection and transmission coefficients. In the 
case of core-sheet interface we have Zc — —Zg and even 
relatively small oscillations of these quantities result in 



strong oscillations of TZcs f^d 7^^ around the averaged val- 
ues typical for the sheet layer of infinite thickness. Both 
the reflected and transmitted waves at the (cs) interface 
come back to this interface after a finite number of cycles 
returning the acoustic energy which has already been ra- 
diated in both directions. This returned energy results in 
the modified growth rate of the instability and again the 
growth rate influence the reflection and transmission co- 
efficients at the (cs) interface. Because of the multiplicity 
of interfaces, many waves interfere within the core-sheet 
jet configuration and the detailed description of these in- 
teractions seems to be a difficult task. 

We explored different sets of values for the densities, 
Lorentz factor of the core and Mach number of the sheet, 
for antisymmetric or symmetric solutions. In all chosen 
cases, \TZae \ remains smaller than 1, implying that no mag- 
nification of the waves occurs at the external interface. 
Variation of the sheet Mach number slightly shifts the os- 
cillating pattern but does not modify the average growth 
rate as long as the sheet keeps a classical velocity. This 
statement is valid also in the case of sheet velocity exceed- 
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Fig. 4. Same as Fig. 2 and 3 but here Vc = 0.1 and Vs = 0.001 



ing the sheet sound speed (Mg > 1) as long as wc assume 
that the core is relativistic and sheet is nonrelativistic. 
The main source of the instability comes from the inter- 
nal interface where amplification of the waves can take 
place. Generally speaking, the instability develops due to 
the interaction between the core and the sheet, with the 
oscillatory modification superimposed by the presence of 
the external (se) interface. 

For small kg^ , when the effective thickness of the sheet 
is smaller than the half- wavelength, namely Ag > 2{R — 
1)/cosq; the growth rate for the stratified case resembles 
to the solution obtained for a single jet without sheet {R = 
1) or appears as a kind of compromise between the two 
extreme cases of no sheet or infinite sheet. 

For all studied sets of parameters, the temporal growth 
rate of the stratified jet instability tends to the solution ob- 
tained for the case of an infinite sheet or oscillates around 
it for large kg^ (corresponding to large k^). Oscillations 
are always reduced while increasing the sheet thickness. 
Practically all solutions that we obtained converge to the 
case of infinite sheet as soon as R reaches the order of 



5 and often even less. So the presence of a sheet with a 
thickness of only a few times the core component radius 
imposes the growth rate of the instability. It is the physi- 
cal parameters of the sheet relatively to the core medium 
which are decisive for the growth of the instability instead 
of the parameters of the external medium. An illustration 
of this is provided by Fig. 2 and Fig. 3 which show almost 
similar growth rate for the stratified scenario (except for 
the inversion of the maxima and minima) in spite of dif- 
ferent values of the density ratio Vf, = Ps/Pe- 

4. Astrophysical implications 

The effect of the sheet that we study in the 2-dimensional 
slab geometry corresponds to the effect due to an enve- 
lope, a sheath or a cocoon aroimd the jet in the cylindri- 
cal description. The general result that properties of the 
sheet essentially influence the linear growth of the Kelvin- 
Helmholtz instability have several consequences. 

It shows that the maximal growth rate and corre- 
sponding wavenumber can be modified by the presence of 
even a modest envelope. This is illustrated by the case of 
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Fig. 5. The physical interpretation of the effect of oscillating 
growth rate is as follows: The sound wave emitted at the in- 
ternal interface first reflects at the external interface and then 
returns to the internal interface and interferes with the new 
emitted wave. The waves drawn within the sheet layer demon- 
strate the integer number of wave periods between interfaces, 
not the way of deformation of sheet (the acoustic waves are 
longitudinal waves). The case L = A3 is shown. 



FM, R = 5 
0.10 1 




0.0 0.2 0.4 0.6 0.8 1.0 



Fig. 6. The maxima of the growth rate coincide with the inte- 
ger multiples of wavelength on the path of the acoustic wave. 
We show the growth rate versus wavenumber k in the linear 
scale for the case presented in Fig. 4 (fundamental mode, R 
= 5). The high bars located below the curve of growth rate 
indicate the fulfillment of the resonance condition (p3) and the 
small bars indicate the fulfillment of the antiresonance condi- 
tion (H) 



Fig. 4, where an envelope not very different from the core- 
component (same thickness with R = 2, density ratio of 
Pc/ps = 0.1) induces an increase of w™'^'^ from 0.03 to 0.2 
and of k^'^^ from 0.1 to 1 for the fundamental mode. The 
linear growth of the helical mode of such a jet is there- 
fore much faster and the resonant wavelength decreases 
to about 6i?c- Such modification due to the new complex- 
ity driven by the envelope-effect can help to solve some 
difficulties encountered when trying to deduce parameters 
such as densities and temperatures from the radio mor- 
phology of jets, under the assumption that knots and wig- 
gles are respectively due to the growth of pinching and 
helical modes (Ferrari et al, 1983; Zaninetti, Van Horn, 
1988). Indeed, the parameters deduced in such studies for 
the external medium can be considered as applying only 
to an envelope surrounding the jet, with density and tem- 
perature possibly quite different from values deduced for 
instance from X-rays data on intergalactic gas. However 
the model of core-sheet jet is complex and involves much 
more free parameters which makes any model-dependent 
predictions more difficult. 

The determination of physical parameters of astro- 
physical jets based on the model of Kelvin - Helmholtz 
instability is a challenge, but it is also made difficult by 
our lack of detailed knowledge on the nonlinear evolution 
of the instability. Several authors investigated the nonlin- 
ear regime by means of numerical simulations (e.g. Hardee 
et al. 1992, Bodo et al. 1994) in the case of nonrelativistic 
jets and applying analytical methods (Hanasz 1995) in the 
case of relativistic jets. The last approach shows that in 
the phase preceding the formation of shock waves within 
the relativistic jet, the growth of instability slows down 
while the perturbation amplitude grows and finally the 
growth is stopped for the amplitude of lateral displace- 
ment of the jet which is only a fraction of the jet radius. 
The physical reasons responsible for this behaviour can 
be summarized as a strong nonlinear growth of modulus 
of the acoustic impedances of both the internal (relativis- 
tic) and the surrounding (nonrelativistic) gases, associated 
with the growth of perturbation amplitude. The properties 
of the surrounding gas and the presence or absence of the 
sheet layer are important for the properties of the Kelvin- 
Helmholtz instability in this aspect as well. The numerical 
simulations allow to trace the evolution of perturbations 
up to the jet disruption. As Bodo et al. (1994) shown the 
time interval between the formation of first shock waves 
and the beginning of the mixing phase is dependent on 
the values of the Mach number of the jet and the density 
contrast. According to this work the smaller density and 
higher Mach number jets can survive in the form of lam- 
inar flow for a longer time. We would like to point out 
that in analogy to our previous results we can expect that 
the parameters of sheet in the case of stratified jets will 
control the life time of the laminar flow. 

Coming back to the linear approximation we point out 
that the presence of an envelope underdense relatively to 
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the external medium tends to induce faster growth of the 

instability while ovcrdcnsc envelope slows it down. The 
high stability of a large number of radio jets can there- 
fore be better understood if dense and cool envelopes are 
commonly present, as expected in "multilevel" or two- 
component jet models (Smith, Raine, 1985; Baker et al, 
1988; Sol et al, 1989; Achatz et al, 1990). This could ex- 
plain for instance the remarkable stability of the inner jet 
of the wide-angle tail (WAT) radio source 3C465. This 
source is hosted by the central galaxy of the cluster of 
galaxies AbcU 2634. A hot intcrclustcr gas makes the ex- 
ternal medium in which the radio jet propagates. X-ray 
observations provided estimates of its temperature and 
density with distance to the core, Te ~ 2 x lO^'K and 
rie ~ 8 X 10~4(p+e-) cm'^ at 20 kpc (Eilek et al, 1984). 
Prom Faraday rotation measures, Leahy (1984) suggested 
the existence of thermal plasma inside the jet with a den- 
sity of 10~^(p+e^) cm^'^ in the inner hot spots at about 
20 kpc from the core and of 5 x 10~^(p"'"e~) cm~^ on av- 
erage in the tails. We tentatively identify this thermal 
plasma with the sheet component of our study (indeed 
in our description, the slow plasma of the sheet can also 
be present at z < Rc, with the fast-moving plasma of the 
core just streaming through it). In two-flow models, the 
slow component (here the sheet or envelope) is made by a 
collimatcd wind coming from the accretion disc. Assuming 
that its temperature is comparable to or less than the one 
of the gas in the disc at the basis of the outflow, one can 
expect a temperature in the range 10"* — lO'^K for the en- 
velope, from the temperature variation proposed by Phin- 
ncy (1989) in the accretion discs. Taking Tg^lO^ required 
X lG^^(p+e^)cm^''' in the region of the inner hot 
spots to ensure pressiirc equilibrium with the intercluster 
gas at the external interface. Such values of Ug are quite 
in agreement with the Faraday rotation analysis and cor- 
respond to an overdense envelope which tends to stabilize 
the inner jet structure. 

Such type of envelopes could be present in many other 
extragalactic radio sources as proposed in an application 
of two-component jet models (Dole et al, 1995) and as 
suggested by the "tomography" analysis of radio data 
(Rudnick, 1995). Direct estimates of the envelope den- 
sities are not yet available. However the sheath detected 
by Katz-Stone and Rudnick (1994) aroimd the jet of the 
FRII source Cygnus A corresponds to an enhancement of 
the density of radiating particles, and therefore likely to 
an increase of the thermal plasma density as well. Over- 
dense sheaths would influence the dynamics of FRII jets 
and increase their stability. A somewhat analogous role 
was already known to be played by FRII cocoons. Indeed 
numerical simulations emphasized the importance of a co- 
coon for the dynamics of FRII jets (Norman et al, 1984a, b) 
and identified two regimes, mode-dominated or cocoon- 
dominated. Our study deal with laminar flows and can 
not precisely describe FRII-cocoons which are turbulent. 
However it proposes some theoretical interpretation of the 



cocoon-dominated regime observed in numerical experi- 
ments. Our results are analogous to the results obtained 
by Hardee & Norman (1990), who notice that the jet dy- 
namics is determined by properties of the jet and lobe, 
but not by the properties in the undisturbed medium. In 
their numerical simulations the lobe surrounds the jet in a 
manner resembling to our sheet or envelope. In fact, four 
different layers can be expected in such cases, namely the 
inner jet, the sheath, the cocoon and the external me- 
dia. Prom our results on core-sheet jets, we tentatively 
infer that the stability properties of this highly stratifled 
configuration will be dominated by the Kelvin-Helmholtz 
instability of the most unstable interface, likely the inter- 
nal one if the inner jet is relativistic, slightly modified by 
the multiple reflections of the acoustic waves at the other 
interfaces. 

The question of the origin of sheaths or envelopes 
around extragalactic jets is still open. They can just em- 
anate from the inner jets as the result of particle diffusion. 
This could concern the slow-moving radio "cocoon" of the 
jet in 3C273. The interaction itself between the jet and 
its surroundings induces turbulent sheared and entrain- 
ment layers which can lead to the formation of sheaths 
for instance in FRI sources. In FRII sources, the backflow 
from the hot spots constitutes some kind of sheaths as 
well, as already mentioned. Other possibilitic;s proposed 
by two-component models for jets are to generate strat- 
ified jets, inner cores and surrounding sheaths, directly 
from the central engines. In all these cases, the presence 
of an envelope with thickness comparable or slightly larger 
than the central jet radius deeply modifies the interaction 
of the jet with the ambient medium of the host galaxy, 
and somewhat isolates the jet from its surroundings. This 
makes then rather problematic to explain the dichotomy 
between radio-loud and radio-quiet objects by the differ- 
ence in the media where the jets have to propagate. Such 
an argument of difference in the host galaxy interstellar 
medium is often proposed to account for the fact that 
well-developed radio jets are observed in elliptical galaxies 
but not in spiral and Seyfert galaxies, based on the idea 
that jets are rapidly destroyed by their interaction with 
the dense and inhomogeneous interstellar gas expected in 
spiral galaxies. In fact, the interstellar medium can in- 
fluence the radio source evolution in two ways, flrst by 
its direct interaction with the jet, second by its primor- 
dial role during the accretion process. Our results favour 
the latter option since the former one appears easily in- 
hibited by specific internal jet properties. Within such a 
view, the dichotomy between radio-loud and radio-quiet 
objects should arise from specific properties of the central 
engine possibly acquired during its formation and accre- 
tion phase. 
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